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Abstract 
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are obtained by means of the Nikiforov-Uvarov method, in closed form. This physical choice of the 
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I. INTRODUCTION 



Within the framework of the Dirac equation the spin symmetry arises if the magnitude of 
the attractive scalar potential S{r) and repulsive vector potential are nearly equal, S{r) ~ 
V{r) in nuclei {i.e., when the difference potential A(r) = V{r) — S{r) — A — constant). 
However, the pseudospin symmetry occurs when S{r) ~ —V{r) are nearly equal {i.e., when 
the sum potential E(r) = V{r) + S{r) = A = constant) [1-3]. The spin symmetry is relevant 
for mesons [4]. The pseudospin symmetry concept has been applied to many systems in 
nuclear physics and related areas [2-7] . It has also been used to explain features of deformed 
nuclei [8] , the super-deformation [9] and to estabhsh an effective nuclear shell- model scheme 
[5,6,10]. The pseudospin symmetry introduced in nuclear theory refers to a quasi-degeneracy 
of the single-nucleon doublets and can be characterized with the non-relativistic quantum 
numbers (n, l,j = l + 1/2) and {n — 1,1 + 2, j = I + 3/2), where n, I and j are the single- 
nucleon radial, orbital and total angular momentum quantum numbers for a single particle, 
respectively [5,6]. The total angular momentum is given as j — I + s, where I — I + 1 is a 
pseudo- angular momentum and 's — 1/2 is a pseudospin angular momentum. In real nuclei, 
the pseudospin symmetry is only an approximation and the quality of approximation depends 
on the pseudo-centrifugal potential and pseudospin orbital potential [11]. In ref. [12], 
Alhaidari et al. investigated in detail the physical interpretation on the three-dimensional 
Dirac equation in the presence of exact spin symmetry limitation A(r) = and pseudospin 
symmetry limitation E(r) = 0. 

Some authors have applied the pseudospin symmetry on several physical potentials, such 
as the harmonic oscillator [12-16], the Woods-Saxon potential [17], the Morse potential [18- 
20], the Hulthen potential [21], the Eckart potential [22-24], the molecular diatomic three- 
parameter potential [25], the Poschl- Teller potential [26] and the Rosen-Morse potential 
[27]. 

On the other hand, the problem of the spatially-dependent effective mass is presenting a 
growing interest along the last few years [28-31]. Many authors have used different methods 
to study the partially exactly solvable and exactly solvable Schrodinger, Klein-Gordon (KG) 
and Dirac equations in the presence of variable mass having a suitable mass distribution 

functions in ID, 3D and/or any arbitrary D-dimensional cases for various type potentials 
[32-45]. In the context of spatially-dependent mass, we have used and applied a recently 
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proposed approximation scheme [45] for the centrifugal term to find a quasi-exact analytic 
bound-state solution of the radial KG equation with spatially-dependent effective mass for 
scalar and vector Hulthen potentials in any arbitrary dimensional space D and orbital an- 
gular momentum quantum number / within the framework of the Nikiforov-Uvarov (NU) 
method [46,47]. In addition, the s-wave bound state solution of the (H-l)-dimensional KG 
equation with mass inversely proportional to the distance from the force center for the in- 
versely linear potential was obtained [48]. Two particular cases were studied, the exact 
spin symmetry when the vector potential and the scalar potential are equal in magnitude 
S{r) = V{r) and the pseudospin symmetry when the vector potential is equal to the scalar 
potential in magnitude but not in sign S{r) = —V{r). Very recently, the study of rcf. [48] has 
also been extended by us [49] to study the bound-state solutions of the (3 + l)-dimensional 
KG equation with position-dependent bosonic mass function m(r) — mo (1 -l- Xobr~^) , where 
r 7^ for the attractive scalar Coulomb-like potential S{r) — —Ksr~^ with Kg — HcQs being 
the coupling constant, taking into consideration the general mixings of scalar and vector 
Lorentz structure potential. 

In the present paper, our aim is to study the exact solution of the spatially dependent 
Dirac equation with the Coulomb-like field potential for any arbitrary spin-orbit quantum 
number n. Under the conditions of the spin symmetry S{r) ~ V{r) and pseudospin symmetry 
S{r) ~ —V{r), we investigate the bound state energy eigenvalues and corresponding upper 
and lower spinor wave functions in the framework of the NU method. We also show that 
the spin and pseudo-spin symmetry Dirac solutions can be reduced to the S{r) — V{r) and 
S{r) = —V{r) Klein-Gordon solutions [49] in the cases of exact spin symmetry limitation 
A(r) = and pseudospin symmetry limitation S(r) = 0, respectively. Further, the solutions 
of the Dirac equation with constant mass can be also generated from the general solution 
when the constant b in the mass function is set to zero. 

The paper is organized as follows. In sect. 2, we outline the NU method. Section 3 is 
devoted for the analytic bound state solutions of the (3-1-1) -dimensional Dirac equation with 
spatially dependent mass function for the quantum system obtained by means of the NU 
method. The spin symmetry and pseudo-spin symmetry solutions arc investigated. In sect. 
3, we study the cases k = ±1 {1 = 1 = 0, s-wave), the constant mass and the non-relativistic 
limit and compare with other wave equations and models. Finally, the relevant conclusions 
are given in sect. 4. 
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II. NU METHOD 



The NU method is briefly outhned here and details can be found in ref . [46] . This method 
was proposed to solve the second-order differential equation of hypergeometric-type: 

m + + ^^Mr) - 0, (1) 

where a{r) and a(r) are polynomials, at most, of second-degree, and r(r) is a first-degree 
polynomial. In order to find a particular solution for eq. (1), let us decompose the wave 
function ■0n('") as follows: 

Mr)^<P{r)yn{r), (2) 

and use 

[(7(r)p(r)]' = r(r)p(r), (3) 

to reduce eq. (1) to the form 

a{r)y';{r) + r{r)y'^{r) + Xy^{r) = 0, (4) 

with 

r(r) = f(r) + 27r(r), /(r) < 0, (5) 

where the prime denotes the differentiation with respect to r. One is looking for a family of 
solutions corresponding to 

A = A„ = -nr'(r)-^n(n-l)a"(r), n = 0,l,2,--., (6) 

The yn(r) can be expressed in terms of the Rodrigues relation: 

where Bn is the normalization constant and the weight function p(r) is the solution of the 
differential equation (3). The other part of the wave function (2) must satisfy the following 
logarithmic equation 

(t)'{r) 7r(r) 

{°) 



4>{r) a{r) 
By defining 

A; = A-7r'(r). (9) 

one obtains the polynomial 



<r) = \ [a'{r) - f (r)] ± ^\ [a'(r) - f(r)]' - 5(r) + ka{T), (10) 

where 7r(r) is a parameter at most of order 1. The expression under the square root sign in 
the above equation can be arranged as a polynomial of second order where its discriminant 
is zero. Hence, an equation for k is being obtained. After solving such an equation, the k 
values are determined through the NU method. 

In this context, we may also derive a parametric generalization from the NU method valid 
for most potentials under consideration. To do this, we begin by writting the hypergeometric 
equation in general parametric form as 

[r (c3 - Cir)f ip'^{r) + [r (cs - dr) (ci - C2r)] ^^(r) + (-^ir^ + - 6) ^n(r) = 0, (11) 

with 

f (r) = ci - C2r, (12) 

cr(r) = r (c3 - C4r) , (13) 

?(r) = -6r' + 6r-e3, (14) 

where the coefficients q {i — 1,2,3,4) and the analytic expressions {j — 1,2,3) are 
calculated for the potential model under consideration. Overmore, comparing eq. (11) with 
it's counterpart eq. (1), we obtain the analytic polynomials, energy equation and wave 
functions together with the associated coefficients expressed in general parameteric form in 
Appendix A. 



III. ANALYTIC SOLUTION OF THE DIRAC-COULOMB-LIKE PROBLEM 

In spherical coordinates, the spatially-dependent mass Dirac equation for fermionic 
massive spin-1/2 particles interacting with arbitrary scalar potential S{r) and the time- 
component V{r) of a four-vector potential can be expressed as [27,50-53] 

[ca-p + 13 {m(ry + S(r)) + V(r) - E] V'„.(r) = 0, Vn«(r) = Vn«(r, 6, 0), (15) 

where E is the relativistic energy of the system, m(r) is the spatially-dependent mass of 
the fermionic particle, p = — iV is the momentum operator, and a and (3 are 4x4 Dirac 
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matrices, which have the following forms, respectively [50-53] 

(16) 




where I denotes the 2x2 identity matrix and a are three-vector Pauh spin matrices 

ai= I ° I , ^2= I ° "M , ^3= I ^ ° I . (17) 





For a particle in a spherical (central) field, the total angular momentum operator J and the 
spin-orbit matrix operator K = (cr • L -|- I) commute with the Dirac Hamiltonian, where 
L is the orbital angular momentum operator. For a given total angular momentum j, the 
eigenvalues of K are k — —{j -\- 1/2) lor aligned spin (si/2, P3/2, ete.) and k — j -\-l/2 lor 
unaligned spin (pi/2, (^3/2, etc.). The spinor wave functions can be classified according to the 
radial quantum number n and the spin-orbit quantum number k and can be written using 
the Pauli-Dirac representation: 

V'n«(r) = - \ , (18) 



r 



lGn.{T)YUe,4>) 



where Fn^ir) and Gn^{r) are the radial wave functions of the upper- and lower-spinor com- 
ponents, respectively, 1^^^(^,0) and Yj^{9,(j)) are the spherical harmonic functions coupled 
to the total angular momentum j and it's projection m on the z axis. The orbital and 
pseudo-orbital angular momentum quantum numbers for spin symmetry / and pseudospin 
symmetry / refer to the upper- and lower-components, respectively. For a given spin-orbit 
quantum number k — ±1, ±2, • • • , the orbital angular momentum and pseudo-orbital an- 
gular momentum are given hy I — \k, + 1/2| — 1/2 and I — \k, — 1/2| — 1/2, respectively. 
The quasi-degenerate doublet structure can be expressed in terms of a pseudo-spin angular 
momentum s = 1/2 and pseudo-orbital angular momentum I which is defined as I = I + 1 
for aligned spin j = 1 — 1/2 and / = / — 1 for unaligned spin j = 1 + 1/2. For example, 
(3si/2, 2^3/2) and (8^1/2,2^3/2) can be considered as pseudospin doublets. 

Substituting eq. (18) into eq. (15), we obtain two radial coupled Dirac equations for the 
spinor components 

ijr + 7) ^'^'^^^^ " + ^"'^ " ^^^^^ ^^^^^ 



G„«(r) = (m(r)c2 - E^, + E(r)) F,,(r), 



(19b) 



d K 
dr r ^ 

where A(r) = V{r) — S{r) and E(r) = V^(r) + S{r) are the difference and sum potentials, 
respectively. Ehminating GnK{r) in eq. (19a) and F„^(r) into eq. (19b), we get two second- 
order differential equations for the upper and lower spinor components as 



K (/€ + 1) 1 

dr"^ h^c^ 
d^ k{k-1) 1 



U-{r)U+{r) + y ^ ^ ' 



U-(r) \dr r 



F„«(r) = 0, (20a) 



Gn.{r) = 0, (20b) 



where U-{r) — m (r) + — A (r) and U+{r) — m(r)c^ — £'„K + S(r), are the difference and 
the sum functions, respectively. Also, y_(r) = - ^ and g+{r) = c^*^ + ^ 

being the derivative of the mass function minus the difference potential and the derivative 
of the mass function plus the sum potential, respectively. From the above equations, the 
energy eigenvalues depend on the quantum numbers n and k, and also the pseudo-orbital 
angular quantum number / according to «;(«; — 1) = /(Z-|-1), which implies that j = / ± 1/2 
are degenerate for i 0. 

At this stage, we take the vector potential in the form of an attractive Coulomb-like field 

as 

E(r) = V{r) = = g, r ^ 0, (21) 

where is being a vector dimensionless real parameter coupling constant and he is being 
a constant with J.fm dimension. Equations (20a) and (20b) can not be solved analytically 
because of the last term in the equations, we find it convenient to solve the mathematical 
relation c^dm{r)/dr = —dE{r)/dr = —dV{r)/dr for the sake of eliminating this term. We 
find out that mass function should be taken as 

m(r) = mo -l- mi/r, (22a) 

mi = moAo&, Aq = h/rrioc, (22b) 

with mo and mi are the integration constant (rest mass of the fermionic particle) and the 
perturbed mass, respectively. Furthermore, b is the dimensionless real constant to be set 
to zero for the constant mass case (i.e., mi = 0) and Ao is the compton-like wavelength in 
fm units. It is worth mentioning that the above choice of the mass function of Coulombic 
form [48,49] is mostly suitable for modeling the well-known pseudo-Coulomb (Kratzer-type) 



potential [54-56] . The interaction field has much impact on the choice of the mass function 
which, in the present case, is inveresely proportional to the distance between the two nuclei 
at short distances m(r) ~ ^ and constant at long distances m(r — >■ oo) ~ mo- 



A. Spin symmetric solution of the Coulomb potential 

In the case of exact spin symmetry S{r) ~ V{r) {dA{r)/dr — 0, i.e., A(r) — A — 
constant), eq. (20a) can be approximately written as 

- - ^ + - A] [m{ry - E^^ + E(r)] } F„.(r) = 0, (23) 

where k, = I for k, < and k = — {I + 1) for k, > 0. The spin symmetric energy eigenvalues 
depend on n and k, i.e., En^ = E{n, k{k+ 1)). In the last equation, the choice of S(r) = 
2V{r) — )• V{r) as mentioned in ref. [12] enables one to reduce the resulting relativistic 
solutions into their non-relativistic limit under appropriate transformations. However, if we 
set E(r) = 2y(r), this allows us to compare the resulting Dirac's solutions with those ones 
derived before in refs. [47-49] regarding the solution of the KG equation for equal mixings 
Sir) = V{r) (i.e., A = 0). 

Substituting eqs. (21) and (22) into eq. (23), allows us to decompose the Dirac equation 
and then leading to obtain the Schrodinger-like equation in the spherical coordinates for the 
upper-spinor component F„i^{r), 



Fn.{r) = 0, (24) 



where we have defined 

en. = ^mlc^ - El^ - Aim^c^ - E^:) > 0, (25a) 

;g = [2g (moc^ + E^,-A)+b{A- 2moc')] , (25b) 

-f^b{b-2q) + K{K + l). (25c) 

For instance, when A = 0, the following constraint Eni < itlqc^ must be fulfilled for bound 
state solutions. The quantum condition is obtained from the finiteness of the solution at 
infinity and at the origin point (i.e., Fn^i^) — Fn.i^) — 0). In order to solve eq. (24) by 
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means of the NU method, we should compare it with eq. (1). The following values for the 
parameters are found as 

f(r)=0, (7(r)=r, 5(r) = -s^^ + - 7. (26) 

Further, inserting these values into eq. (10), we obtain 

1 1 



7r(r ) = - ± - Vl + 47 + 4(/c-/3)r + 4£2^r2. (27) 



The discriminant of the square root must be set equal to zero, i.e., l+47+4(/c— /3)r+4£^^r^ 
0. Hence, we find the following two constants: 

A;i,2 = ^±£„kV1 + 47- (28) 

In this regard, we can find the possible functions for 7r(r) as 

[ i ± [er^^r + |v/rT4^] for h = f3 + £„,^rT4^, , , 

7r(rj = < (^yj 

[ i ± [e„^r - iyrq^] for h^^- en^^/T+I^. 

According to the NU method, one of the four values of the polynomial 7r(r) is just proper 
to obtain the bound states because r(r) has a negative derivative. Therefore, the selected 
forms of 7r(r) and k take the following particular values 

7r(r) = ^ (1 + \/l + - £nnr, k = (3 - £^^0^4^, (30) 

to obtain 

T(r) = 1 + Vl + 47 - 25n.r, r'(r) = -2£n« < 0, (31) 

where prime denotes the derivative with respect to r. In addition, after using eqs. (6) and 
(9) together with the assignments given in eqs. (26), (30) and (31), the following expressions 
for A are obtained as 

A„ = A = 2n£„«, n = 0, 1, 2, • • • , (32) 
A = ^ - £„K (1 + aA+I^) . (33) 

Now, taking A„ = A, we can solve the above equations to obtain the energy equation for the 
Coulomb-like potential with spin-symmetry in the Dirac theory. 



q {rriQC^ + £;„^ - A) + b - mpc^) 



(34) 
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where 



1 + ^) +b{b-2q)-^. (35) 



For spatially-dependent mass case, i.e., 6 7^ 0, the bound state solutions of the system 
are determined by means of the parameters q and b. It is not difficult to conclude that all 
bound-states appear in pairs, two energy solutions are valid for the particle — and 
the second one corresponds to the anti-particle energy E"- = E~^ in the Coulomb-like field. 

Let us now find the corresponding eigenfunctions for this system. Using eqs. (3) and (8), 
we find 

p{r) = r2'^+^e-2^"''^ (36) 
0(r) = r-^+^e-""'^'', (37) 
and further substituting eqs. (26) and (36) into eq. (7), we get 

yn{r) = A„r-(2^+i)e^--|;^ [,26+n+i^-2.„..] ^ L^^+i(2£„.r), (38) 

where L^{x) is the generalized Laguerre polynomials. The relation Fn^{r) = 4>{r)yn{r) gives 
the radial upper spinor wave function: 

Fn.{r) = Mr'^'e-^-^'Ll'+\2e^,r). (39) 

The above upper-spinor component satisfies the restriction condition for the bound states, 
i.e., 6 > and e„K > 0. Using the normalization condition u{r)^dr = 1 and the or- 



a+l„—x 



thogonality relation of the generalized Laguerre polynomials 
{2n + a + 1) the normahzing factor J\f can be found as [57-59] 



dx — 



^' 2{n + S + l)T{n + 2d + 2y ^ ^ 

where and 5 are given in eqs. (25a) and (35), respectively. 

On the other hand, it is worth noting that if we take the particular values of the coefficients 
displayed in Table 1, then the expressions given in Appendix A can be calculated for eq. 
(24) and tested with the ones calculated before in a very simple way. 

Let us now study some special cases of much concern to the reader. The following cases 
are found to be identical with the solution of the spatially-dependent mass and constant 
mass KG equation for the case spin symmetry case S{r) — V{r) (A(r) = E(r) = 0) [29,51]. 
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(1) If we consider the case when S{r) — V{r), i.e., Qs — Qv — q or A — 0, then eq. (34) 
can be reduced to the following forms: 



and 



UK 



q{b-q) + Bn^^Bl^-h{h-2q) 



q^ + B, 



2 

UK 



q (b-q)- B„iy/Bl^-b{b-2q) 



q^ + B^ 

y ' riK 



(41a) 



(41b) 



where is defined by 



B„ 



n 



k) +b{b-2q)\ . 



(42) 



Obviously, the bound-state solutions of the particle and anti-particle do exist. 
(2) When the mass is constant, i.e., 6 = (mi = 0), we have 



EL = i"^"^!!^ . 'I rnoc^ and - ^ 



(n + K + l)^ + g2 
which gives 



mo<? and — — moc^, n = 0, 1, 2, ■ ■ ■ , k — ±1, ±2, 



(43) 



UK 



(n+/+l)2 g2 UK 

T—j-wT—omQC , K < u 



-moc^, / = 0, 1, 2, 



(44) 



where n and k signify the usual radial and spin-orbit quantum numbers. The particle has 
bound state solution whereas anti-particle has continuum solution for all states. Further, 
the wave functions for particle turn to be 



X 



L2«+i (^1-^^ mlc^ - Ef^T^ , K>0, EL< moc^ (45) 

whereas the anti-particle has no wave functions. 

(3) When the potential coupling constant is set to g = 6/2, the spectrum of the spatially- 
dependent mass Dirac particle in potential field qs — qv — q is similar to the spectrum of 
constant mass Dirac particle in the potential fields qs — —qv, that is. 



(46a) 
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{n + K + lf + 

(4) When the potential couphng constant is set to = 6/2, the spectrum of the spatially- 
dependent mass Dirac particle in potential field Qg — — g„ is similar to the spectrum of 
constant mass Dirac particle in the potential fields Qs = Qv = Q, that is, 

{n + K + 1) + 

E^i = -moc'. (47b) 



B. Pseudospin symmetric solution of the Coulomb potential 



In the presence of the pseudospin symmetry S{r) ~ —V{r) {i.e., dT,{r)/dr = 0, or 
S(r) — A — constant), eq. (20b) can be exactly written as 



(f _ K (k - 1) 



^ [m{ry + E^, - A(r)] [m{ry - E^, + A] | C(r) = 0, (48) 



where the energy eigenvalues En^ depend only on n and k,, i.e., E^^ = E{n, k{k — 1)). We 
take the difference of the potential as 2V{r) which allows us to compare our final results 
with the ones calculated for the KG equation derived in ref. [49]. Substituting 



A(r) = V{r) 



hcQy 



together with eq. (22) into the last equation, we obtain a Schrodinger euation for the lower 
component Gnnij'), 



dr"^ 



Gnn{r) = 0, GnM = Gn.{^) = 0, 



(49) 



where we have defined 



SuK = ^ mlc^ -El^ + A (moc2 + K«) > 0, 



(50a) 



^ = -— [2q (moc" -E„^ + A)+b (2moc' + A)] , 



(50b) 

7 = 6 (6 + 2g) + K (k + 1) . (50c) 

To avoid repetition in the solution of eq. (49), a first inspection for the relationship be- 
tween the present set of parameters (£„k, (3, 7) and the previous set {snK, 7) tells us that 
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the negative energy solution for pseudospin symmetry, where S{r) — —V{r), can be ob- 
tained directly from those of the positive energy solution above for spin symmetry using the 
parameter map [50-52]: 

F„«(r) O Gnn{r),V{r) -V{r) (or q^-q), -E^,^ and A -A. (51) 

Following the previous results with the above transformations, we finally arrive at the energy 
equation, 

2 

, (52) 



mlc" -El^ + A (moc^ + E^^) 



q (moc^ - Enn + -4) + h {nior^ + |) 



n+l + J{n-^y + b{b + 2q) 



and lower-spinor component wave function. 



2»7+3 



2{n + 7] + l)T{n + 2rj + 2) 



r^+^e-^-^'Ll^+\2e,,,r), 



(53) 



where V — y {^^ ~ ^) + b(b + 2q) — ^. The finiteness of the above wave functions for the 
bound states is achieved once 77 > and > 0. Let us investigate the following specific 
cases of interest. 

(1) Considering the case when S{r) — —V{r), i.e., qs — —qv — q or A — 0, then eq. (52) 
can be reduced to the following forms: 



-q {b + q) + Bn.^Bl^-b{b + 2q) 



-q (b + q)- Br.i^Bl,-b{b + 2q) 



moc 



q^ + B^ 



2 

UK 



moc 



(54a) 



(54b) 



where Bn^ is defined by 



B„ 



+ +b{b + 2q) I . 



(55) 



We note that the bound state solutions of the particle and anti-particle are available. 
(2) When the mass is constant, i.e., 6 = 0, then we obtain 



riK 



{n + Kf-q^_ 2 



(n + k) + q^ 



moc, El^ = -moc, n = 0, 1, 2, • • • , k = ±1, ±2, 



(56) 
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and the lower-spinor wave functions: 



2k+1 



-r e 



where k = 1, 2, 3, • • ■ , and E^k < itiqc^. 



IV. DISCUSSIONS 



(57) 



In this section, we are going to study two special cases of the energy eigenvalues given 
by eqs. (34) and (52). First, let us study s-wave case I — Q {k — —1) and I — Q {k — 1) case 

2 



mlc^ - El_^ - A (moc' - = 



m^c^ - + A (moc' + K,i) 



2q {moc^ + En-i - A) +b{A- 2moc^ 
2n + 1 + + 46(6 - 2g) 



2q {moc^ - En,i + A) + b{A + 2moc^ 



(58) 



(59) 



2n + 1 + ^/TTWbT2q) 
For constant mass case, if one set A = and 6 = into eqs. (58) and (59), we obtain for 
spin and pseudo-spin symmetric Dirac theory. 



En,-1 — 



EnA — 



r 



-rriQC , 



{n+lf-q^_ , 



-rriQC . 



(60a) 



(61a) 



(n + l)' + g2 

respectively, which are identical with the s-wave results in ref. [48] for S{r) = V{r) and 
S{r) — —V{r), respectively. Let us now discuss the non-rclativistic limit of the energy 
eigenvalues and wave functions of our solution for the spatially dependent mass and for 
the constant mass Schrodinger equation. If we take A — and put 2g — > g [12], i.e., 
S{r) — V{r) — S(r), the non-relativistic limits of energy equation (34) and wave functions 
(39) under the following transformations EnK + moC^ ^ 2/x//i^ and E^k — rriQC^ E^i [27,50] 
become 

(g - hf 



E., 



nl 



2^2 



{n-r\ + ^J{\ + lY + h{h-2q)^ 



2 ' 



(62) 
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Fniir. 
respectively, where 



.ir , 



(63) 



y'(2;+l)2+46(6-q)+2 



^ (^2n + l + ^(2Z + l)' + 46(6-g)j T + 1 + ^ {21 + if + Ab{b - q) 



(64) 



Furthermore, when mass becomes constant (6 = 0), we find the well-known non-relativistic 
energy spectrum and wave functions of the Coulomb problem, 

„2 



2/12 (n + Z + 1) 



(65) 



which are identical to those given in refs. [60-62]. 



n\ (2^-fK, 



2i+3 



2(n + / + l)r(n + 2/ + 2) 

(66) 



V. CONCLUSIONS 



To summarize, we have presented the bound state solutions of the spatially-dependent 
mass Dirac equation with the Coulomb-like field potential under the conditions of the spin 
symmetry and pseudospin symmetry. We have obtained an explicit expressions for energy 
eigenvalues and associated wave functions for arbitrary spin-orbit k state with a specific 
choice for the mass function that provides an exact solution to the pseudospin symmetric 
Dirac equation and an approximate solution for the spin symmetric Dirac equation. This 
suitable choice of mass function which is in the form of Coulomb-like potential enables us 
to solve the spatially-dependent Dirac equation analytically and reduces it to the constant 
mass solution as well as when 6 = (mi = 0). The most stringent interesting result is 
that the present spin and pseudo-spin symmetries can be easily reduced to the previously 
found Klein-Gordon solution once S{r) — V{r) and S{r) — —V{r) {i.e., when A — 0), 
respectively. The resulting solutions of the wave functions are being expressed in terms of 
the generalized Laguerre polynomials. Obviously, when the coupling potential parameters 
are adjusted to some specific values, particularity when q — 6/2, the spectra of the mass 
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varying Dirac particle for the case Qs — Qv {Qs — —Qv) turn to become similar to the spectra 
of the constant mass Dirac particle for the case Qs — —Qv {Qs — Qv), respectively. In the 
limit of constant mass (6 = 0), the solution for the energy eigenvalues and wave functions 
are reduced to those ones given in literature. Also, when spin-orbit quantum number k = 0, 
the problem reduces to s- waves solution as in [48]. 
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Appendix A: Parametric Generalization of the NU Method 

Our systematical derivations hold for any potential form, 
(i) The analytic results for the key polynomials [47] : 

7r(r) = C5 + ^ - — {ciyfc^ + ^/c^ - C3C6) r, (Al) 
C3 

-\ (C3C8 + 2C4C9 + 2VC9C10) , (A2) 
2 

T(r) = C3 + 2^ (C3C4 + + (A3) 

ca 

2 

r'(r) = (C3C4 + C4v^ + v^) < 0. (A4) 

C3 

(ii) The energy equation: 

C2n - (2n + 1) C6 + — (2n + 1) {^/c^ + C4v/c^) 
C3 

+ n (n - 1) C4 + ^ (C3C8 + 2C4C9 + 21/C9C10) = 0. (A5) 
C3 

(iii) The wave functions: 

p(r) =r'^"(c3-C4r)^i% (A6) 

0(r) = r^i=^(c3 - C4r)^", (A7) 

y„(r) = Pi'^"''^-)(c3-C4r), (A8) 

V^n(r) = (t){r)yn{r) = Ar'"(c3 - C4r)'="Pi^"''=")(c3 - 2c4r), (A9) 

where Pn'^\c^ — C/^r) are the Jacobi polynomials and An is a normalizing factor. 

When C4 = 0, the Jacobi polynomial turn to be the generalized Laguerre polynomial and 

the constants relevant to this polynomial change are 

lim P^'^"'^i^)(c3 - 2c4r) = ^^"(cisr), (AlO) 

C4— 

lim (C3 - C4r)^" = e'^^''^ (All) 

where L^"(ci5r) are the generalized Laguerre polynomials and Bn is a normalizing factor, 
(iv) The relevant coefficients q (i = 5, 6, • • • , 16) are given as follows: 

1 1 

C5 = 2 (^3 - Cl) , C6 = - (C2 - 2C4) , C7 = Cg + {1, (A12) 



17 



Cs = 2C5C6 - ^2, Cg = C5 + ^3, 



Cio — C4 (C3C8 + C4C9) + C3C7, 

2 2 

Cll = — a/c^, C12 = a/cI^' 

Cs C3C4 



(C5 + Vc^) ) Ci4 = (\/ciO - C4C5 - C3C6) , 

C3C4 

_ 2 _ Ci5 

Cl5 — — V^IO) C16 — — . 

C3 2 



(A13) 
(A14) 

(A15) 

(A16) 
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TABLE I: The specific values for the parametric constants necessary for calculating the energy 
eigenvalues and eigenfunctions of the spin symmetry Dirac wave equation (23). 



Constant Analytic value Constant Analytic value 
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